Disentangling Crocus Series Verni and Its Polyploids.

Spring crocuses, the eleven species within Crocus series Verni (Iridaceae), consist of di- and tetraploid cytotypes. Among them is a group of polyploids from southeastern Europe with yet-unclear taxonomic affiliation. Crocuses are generally characterized by complex dysploid chromosome number changes, preventing a clear correlation between these numbers and ploidy levels. To reconstruct the evolutionary history of series Verni and particularly its polyploid lineages associated with C. heuffelianus, we used an approach combining phylogenetic analyses of two chloroplast regions, 14 nuclear single-copy genes plus rDNA spacers, genome-wide genotyping-by-sequencing (GBS) data, and morphometry with ploidy estimations through genome size measurements, analysis of genomic heterozygosity frequencies and co-ancestry, and chromosome number counts. Chromosome numbers varied widely in diploids with 2n = 8, 10, 12, 14, 16, and 28 and tetraploid species or cytotypes with 2n = 16, 18, 20, and 22 chromosomes. Crocus longiflorus, the diploid with the highest chromosome number, possesses the smallest genome (2C = 3.21 pg), while the largest diploid genomes are in a range of 2C = 7-8 pg. Tetraploid genomes have 2C values between 10.88 pg and 12.84 pg. Heterozygosity distribution correlates strongly with genome size classes and allows discernment of di- and tetraploid cytotypes. Our phylogenetic analyses showed that polyploids in the C. heuffelianus group are allotetraploids derived from multiple and partly reciprocal crosses involving different genotypes of diploid C. heuffelianus (2n = 10) and C. vernus (2n = 8). Dysploid karyotype changes after polyploidization resulted in the tetraploid cytotypes with 20 and 22 chromosomes. The multi-data approach we used here for series Verni, combining evidence from nuclear and chloroplast phylogenies, genome sizes, chromosome numbers, and genomic heterozygosity for ploidy estimations, provides a way to disentangle the evolution of plant taxa with complex karyotype changes that can be used for the analysis of other groups within Crocus and beyond. Comparing these results with morphometric analysis results in characters that can discern the different taxa currently subsumed under C. heuffelianus.

Ancient Artworks and Crocus Genetics Both Support Saffron's Origin in Early Greece.

Saffron crocus (Crocus sativus) is a male-sterile, triploid flower crop, and source of the spice and colorant saffron. For over three millennia, it was cultivated across the Mediterranean, including ancient Greece, Persia, and other cultures, later spreading all over the world. Despite saffron crocus' early omnipresence, its origin has been the matter of a century-old debate, in terms of area and time as well as parental species contribution. While remnants of the ancient arts, crafts, and texts still provide hints on its origin, modern genetics has the potential to efficiently follow these leads, thus shedding light on new possible lines of descent. In this review, we follow ancient arts and recent genetics to trace the evolutionary origin of saffron crocus. We focus on the place and time of saffron domestication and cultivation, and address its presumed autopolyploid origin involving cytotypes of wild Crocus cartwrightianus. Both ancient arts from Greece, Iran, and Mesopotamia as well as recent cytogenetic and comparative next-generation sequencing approaches point to saffron's emergence and domestication in ancient Greece, showing how both disciplines converge in tracing its origin.

Plant migration under long-lasting hyperaridity - phylogenomics unravels recent biogeographic history in one of the oldest deserts on Earth.

The post-Miocene climatic histories of arid environments have been identified as key drivers of dispersal and diversification. Here, we investigate how climatic history correlates with the historical biogeography of the Atacama Desert genus Cristaria (Malvaceae). We analyze phylogenetic relationships and historical biogeography by using next-generation sequencing (NGS), molecular clock dating, Dispersal Extinction Cladogenesis and Bayesian sampling approaches. We employ a novel way to identify biogeographically meaningful regions as well as a rarely utilized program permitting the use of dozens of ancestral areas. Partial incongruence between the established taxonomy and our phylogenetic data argue for a complex historical biogeography with repeated introgression and incomplete lineage sorting. Cristaria originated in the central southern part of the Atacama Desert, from there the genus colonized other areas from the late Miocene onwards. The more recently diverged lineages appear to have colonized different habitats in the Atacama Desert during pluvial phases of the Pliocene and early Pleistocene. We show that NGS combined with near-comprehensive sampling can provide an unprecedented degree of phylogenetic resolution and help to correlate the historical biogeography of plant communities with cycles of arid and pluvial phases.

Phylogeography of Artemisia frigida (Anthemideae, Asteraceae) based on genotyping-by-sequencing and plastid DNA data: Migration through Beringia.

Artemisia frigida is a temperate grassland species that has the largest natural range among its genus, with occurrences across the temperate grassland biomes of Eurasia and North America. Despite its wide geographic range, we know little about the species' distribution history. Hence, we conducted a phylogeographical study to test the hypothesis that the species' distribution pattern is related to a potential historical migration over the 'Bering land bridge'. We applied two molecular approaches: genotyping-by-sequencing (GBS) and Sanger sequencing of the plastid intergenic spacer region (rpl32 - trnL) to investigate genetic differentiation and relatedness among 21 populations from North America, Middle Asia, Central Asia and the Russian Far East. Furthermore, we identified the ploidy level of individuals based on GBS data. Our results indicate that A. frigida originated in Asia, spread northwards to the Far East and then to North America across the Bering Strait. We found a pronounced genetic structuring between Middle and Central Asian populations with mixed ploidy levels, tetraploids in the Far East, and nearly exclusively diploids in North America except for one individual. According to phylogenetic analysis, two populations of Kazakhstan (KZ2 and KZ3) represent the most likely ancestral diploids that constitute the basally branching lineages, and subsequent polyploidization has occurred on several occasions independently. Mantel tests revealed weak correlations between genetic distance and geographical distance and climatic conditions, which indicates that paleoclimatic fluctuations may have more profoundly influenced A. frigida's spatial genetic structure and distribution than the current environment.

Resolving relationships in an exceedingly young Neotropical orchid lineage using Genotyping-by-sequencing data.

Poor morphological and molecular differentiation in recently diversified lineages is a widespread phenomenon in plants. Phylogenetic relationships within such species complexes are often difficult to resolve because of the low variability in traditional molecular loci. Furthermore, biological phenomena responsible for topological incongruence such as Incomplete Lineage Sorting (ILS) and hybridisation complicate the resolution of phylogenetic relationships among closely related taxa. In this study, we employ a Genotyping-by-sequencing (GBS) approach to disentangle evolutionary relationships within a species complex belonging to the Neotropical orchid genus Cycnoches. This complex includes seven taxa distributed through Central America and the Colombian Chocó, and is nested within a clade estimated to have first diversified in the early Quaternary. Previous phylogenies inferred from few loci failed to provide support for internal relationships within the complex. Our Neighbour-net and coalescent-based analyses inferred from ca. 13,000 GBS loci obtained from 31 individuals belonging to six of the seven traditionally accepted Cycnoches taxa provided a robust phylogeny for this group. The genus Cycnoches includes three main clades that are further supported by morphological traits and geographic distributions. Similarly, a topology reconstructed through maximum likelihood (ML) inference of concatenated GBS loci produced results that are comparable with those reconstructed through coalescence and network-based methods. Our comparative phylogenetic informativeness analyses suggest that the low support evident in the ML phylogeny might be attributed to the abundance of uninformative GBS loci, which can account for up to 50% of the total number of loci recovered. The phylogenomic framework provided here, as well as morphological evidence and geographical patterns, suggest that the six entities previously thought to be different species or subspecies might actually represent only three distinct segregates. We further discuss the limited phylogenetic informativeness found in our GBS approach and its utility to disentangle relationships within recent and rapidly evolving species complexes. Our study is the first to demonstrate the utility of GBS data to reconstruct relationships within young (~2 Ma) Neotropical plant clades, opening new avenues for studies of species complexes that populate the species-rich orchid family.

New microsatellite markers for Xerophyta dasylirioides (Velloziaceae), an endemic species on Malagasy inselbergs.

PREMISE: Microsatellite markers were developed for Xerophyta dasylirioides (Velloziaceae), a species endemic to the Malagasy inselbergs, to explore the impact of its island-like distribution on genetic diversity and gene flow. METHODS AND RESULTS: A total of 7110 perfect microsatellite loci were recovered by shotgun sequencing on an Illumina MiSeq platform. Primer pairs were designed for 40 arbitrarily selected loci. Fifteen primer pairs that generated distinct PCR products were used to genotype 80 individuals of X. dasylirioides from three inselberg populations. All markers were polymorphic, revealing two to 17 alleles in the overall sampling. Levels of observed and expected heterozygosity per locus ranged from zero to 1.000 and from zero to 0.850, respectively. Success rates of cross-amplification in 10 additional species of Xerophyta (X. croatii, X. decaryi, X. isaloensis, X. labatii, X. lewisiae, X. pinifolia, X. retinervis, X. setosa, X. spekei, X. tulearensis) ranged from zero to 70%. CONCLUSIONS: Fifteen newly developed microsatellite markers provide a toolkit for assessing population genetic parameters of X. dasylirioides in its unique island-like habitats.

A pre-Miocene Irano-Turanian cradle: Origin and diversification of the species-rich monocot genus Gagea (Liliaceae).

The Irano-Turanian (IT) floristic region is considered an important center of origin for many taxa. However, there is a lack of studies dealing with typical IT genera that also occur in neighboring areas. The species-rich monocot genus Gagea Salisb. shows a center of diversity in IT region and a distribution in adjacent regions, therefore representing a good study object to investigate spatial and temporal relationships among IT region and its neighboring areas (East Asia, Euro-Siberia, Himalaya, and Mediterranean). We aimed at (a) testing the origin of the genus and of its major lineages in the IT region, (b) reconstructing divergence times, and (c) reconstructing colonization events. To address these problems, sequences of the ribosomal DNA internal transcribed spacer (ITS) region of 418 individuals and chloroplast intergenic spacers sequences (psbA-trnH, trnL-trnF) of 497 individuals, representing 116 species from all sections of the genus and nearly its entire distribution area were analyzed. Divergence times were estimated under a random molecular clock based on nrITS phylogeny, which was the most complete data set regarding the representation of species and distribution areas. Ancestral distribution ranges were estimated for the nrITS data set as well as for a combined data set, revealing that Gagea most likely originated in southwestern Asia. This genus first diversified there starting in the Early Miocene. In the Middle Miocene, Gagea migrated to the Mediterranean and to East Asia, while migration into Euro-Siberia took place in the Late Miocene. During the Pleistocene, the Arctic was colonized and Gagea serotina, the most widespread species, reached North America. The Mediterranean basin was colonized multiple times from southwestern Asia or Euro-Siberia. Most of the currently existing species originated during the last 3 Ma.

Saffron (Crocus sativus) is an autotriploid that evolved in Attica (Greece) from wild Crocus cartwrightianus.

Crocus sativus, the saffron crocus, is the source of saffron, which is made from the dried stigmas of the plant. It is a male-sterile triploid lineage that ever since its origin has been propagated vegetatively. Its mode of evolution and area of origin are matters of long-lasting debates. Here we analyzed chloroplast genomes and genome-wide DNA polymorphisms obtained through genotyping-by-sequencing (GBS) to infer the parent and area of origin of C. sativus. These data were complemented by genome size measurements and analyses of nuclear single-copy genes. We could place 99.3% of saffron GBS alleles in Crocus cartwrightianus, a species occurring in southeastern mainland Greece and on Aegean islands, identifying it as the sole progenitor of the saffron crocus. Phylogenetic and population assignment analyses together with chloroplast polymorphisms indicated the C. cartwrightianus population in the vicinity of Athens as most similar to C. sativus. We conclude that the crop is an autotriploid that evolved in Attica by combining two different genotypes of C. cartwrightianus. Triploid sterility and vegetative propagation prevented afterwards segregation of the favorable traits of saffron, resulting in worldwide cultivation of a unique clonal lineage.

Phylogeny of the saffron-crocus species group, Crocus series Crocus (Iridaceae).

Phylogenetic relationships among the taxa of Crocus series Crocus are still unclear, preventing the understanding of species diversity and the evolution of the important spice saffron (Crocus sativus). Therefore, we analyzed sequences of two chloroplast (trnL-trnF, matK-trnK) and three nuclear (TOPO6, ribosomal DNA ETS and ITS) marker regions to infer phylogenetic relationships among all species belonging to series Crocus. Our phylogenetic analyses resolved the relationships among all taxa of the series. Crocus hadriaticus and the former C. pallasii subspecies appeared polyphyletic. The latter deserve elevating the subspecies to species rank, while for C. hadriaticus a detailed study of species boundaries is necessary. Multi-locus and also genome-wide single nucleotide polymorphism data obtained through genotyping-by-sequencing placed C. sativus within C. cartwrightianus with no indication that other Crocus species contributed to the evolution of the triploid. Our analyses thus made an autotriploid origin of C. sativus from C. cartwrightianus very likely.

Development of microsatellite markers and assembly of the plastid genome in Cistanthe longiscapa (Montiaceae) based on low-coverage whole genome sequencing.

Cistanthe longiscapa is an endemic annual herb and characteristic element of the Chilean Atacama Desert. Principal threats are the destruction of its seed deposits by human activities and reduced germination rates due to the decreasing occurrence of precipitation events. To enable population genetic and phylogeographic analyses in this species we performed paired-end shotgun sequencing (2x100 bp) of genomic DNA on the Illumina HiSeq platform and identified microsatellite (SSR) loci in the resulting sequences. From 29 million quality-filtered read pairs we obtained 549,174 contigs (average length 614 bp; N50 = 904). Searching for SSRs revealed 10,336 loci with microsatellite motifs. Initially, we designed primers for 96 loci, which were tested for PCR amplification on three C. longiscapa individuals. Successfully amplifying loci were further tested on eight individuals to screen for length variation in the resulting amplicons, and the alleles were exemplarily sequenced to infer the basis for the observed length variation. Finally we arrived at 26 validated SSR loci for population studies in C. longiscapa, which resulted in 146 bi-allelic SSR markers in our test sample of eight individuals. The genomic sequences were also used to assemble the plastid genome of C. longiscapa, which provides an additional set of maternally inherited genetic markers.

Phylogeny of Crocus (Iridaceae) based on one chloroplast and two nuclear loci: ancient hybridization and chromosome number evolution.

Crocus consists of about 100 species distributed from western Europe and northern Africa to western China, with the center of diversity on the Balkan Peninsula and in Asia Minor. Our study focuses on clarifying phylogenetic relationships and chromosome number evolution within the genus using sequences of the chloroplast trnL-F region, the nuclear ribosomal DNA internal transcribed spacer (ITS) region, and a part of the nuclear single-copy gene pCOSAt103. In a combined dataset of ITS and trnL-F sequences, 115 individuals representing 110 taxa from both subgenera and all sections and series of Crocus were analyzed with Bayesian phylogenetic inference. For pCOSAt103 79 individuals representing 74 Crocus taxa were included, and for the majority of them PCR amplicons were cloned and up to eight clones per individual were sequenced to detect allopolyploidization events. Romulea species were included as outgroup in both analyses. Characteristics of seed surface structures were evaluated by scanning electron microscopy. Phylogenetic analysis of ITS/trnL-F data resulted in a monophyletic genus Crocus, probably monophyletic sections Crocus and Nudiscapus, and inferred monophyly for eight of the 15 series of the genus. The C. biflorus aggregate, thought to be consisting of closely related subspecies, was found to be polyphyletic, the taxa occurring within three major clades in the phylogenetic tree. Cloning of pCOSAt103 resulted in the detection of homoeologous copies in about one third of the taxa of section Nudiscapus, indicating an allotetraploid origin of this section. Reconstruction of chromosome number evolution along the phylogenetic tree using a probabilistic and a parsimony approach arrived at partly contradictory results. Both analyses agreed however on the occurrence of multiple polyploidization and dysploidy events. B chromosomes evolved at least five times independently within the genus, preferentially in clades characterized by karyotype changes.